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BACKGROUND 


The  U.S.  Army  combined  effects  aluminized  explosives  (PAX-29,  PAX-30,  and  PAX-42) 
were  previously  demonstrated  to  achieve  both  excellent  metal  pushing  and  high  blast  energies 
in  both  cylinder  tests  and  warhead  configurations.  The  excellent  metal  pushing  capability  of 
these  combined  effects  explosives  is  due  to  the  earlier  exothermic  conversion  of  aluminum  to 
aluminum  oxide  as  compared  to  traditional  blast  explosives.  However,  the  traditional  Chapman- 
Jouguet  detonation  theory  does  not  explain  the  observed  detonation  states  achieved  by  these 
combined  effects  explosives.  The  eigenvalue  detonation  theory  explains  the  observed  behavior 
and  was  incorporated  into  the  computational  design  process  for  warhead  configurations  (refs.  1 
through  3).  An  analytic  cylinder  test  model  has  long  been  used  by  the  U.S.  Army  Armament 
Research,  Development  and  Engineering  Center  (ARDEC),  Picatinny  Arsenal,  New  Jersey  for 
explosive  equation  of  state  calibration  and  verification.  This  analytic  model  has  been  shown  to 
provide  close  agreement  to  high  rate  continuum  modeling.  However,  the  analytic  model  was 
based  on  adiabatic  expansion  along  the  principle  isentrope  from  the  Chapman-Jouguet  state. 

As  it  has  been  hypothesized  that  the  current  combined  effects  explosives  expand  from  a  weak 
point  on  the  fully  reacted  Hugoniot,  the  validity  of  the  analytic  cylinder  test  comes  into  question 
for  these  explosives.  The  analytic  cylinder  test  model  was  recently  updated  to  include  the 
eigenvalue  detonation  theory  and  associated  adiabatic  expansion  from  the  fully  reacted 
Hugoniot  weak  point. 


TAYLOR  LONG  BOMB  MODEL 


Gurney  formulation  is  a  traditional  analytic  method  used  for  high  explosive  material 
acceleration  modeling  (ref.  4).  The  work  of  Taylor  (ref.  5)  provides  the  framework  for  a  more 
fundamental  methodology  for  modeling  exploding  cylinders,  including  axial  flow  effects.  Based 
on  the  assumption  of  heavy  confinement,  Taylor  proposed  the  use  of  Reynolds  hydraulic 
treatment  and  isentropic  detonation  products  flow  to  model  a  heavy  long  cylindrical  bomb 
expansion.  The  model  is  very  similar  to  the  Delaval  nozzle  theory  (ref.  6),  but  includes  the 
momentum  of  the  cylinder  wall.  A  diagram  of  a  Taylor  heavy  long  cylinder  expansion  due  to  high 
explosive  detonation  is  presented  in  figure  1.  It  should  be  noted  that  flow  velocities  are  relative 
to  the  detonation  velocity,  D.  Using  the  Jones-Wilkens-Lee-Baker  [JWLB  (ref.  7)]  thermo¬ 
dynamic  equation  of  state,  the  Taylor  heavy  long  cylindrical  bomb  model  can  be  represented  as 


follows. 


Mass: 


p0Dr02  =  pUr2 


(1) 


Linear  Momentum: 


2nrP  =  mD  — 

dt 


(2) 


Energy: 


(3) 


1 


Principle  isentrope: 


P 


RjPo 

T.iAie  p  + 


C 


Poy((0+l) 


de  = 


Taylor  angle: 

v  =  2D  sin  — 

2 


(4) 

(5) 

(6) 


Based  on  the  assumption  of  heavy  confinement,  this  model  neglects  radial  flow  effects  in 
the  detonation  products.  Additionally,  the  detonation  products  flow  is  assumed  to  be  isentropic 
from  the  Chapman-Jouguet  state  and  cylinder  strength  is  neglected.  To  achieve  an  easier  com¬ 
putational  form,  the  following  reduction  is  made. 


c  =  p07rr0" 


(1)  =>  (f)2  = 

Vr0/  p 


vr  =  ^  =  DsinQ,  (2)  =>  2nrP 


.mD2dcos0  =  mD225!i^ 


dr 


dr 


(7) 

(8) 

(9) 


2 


(6) 


^L  =  ^LP^d-  =  --d(-)2 

dr  M  2  m  p0  \r0J 


~«W2}-m27:dp 

®3dT  =  icK^)-^‘">] 


(3)  =>^-  +  e0  =-^  +  -±  +  ecj  =  V-  +  P-  +  e 

2  u  2  pCj  CJ  2  p 


(5)  =>  de  =  -Pd  (-) , (13)  =>  d  —  =  -d-  +  Pd-  =  —  -  dP 

W  2  p  p  p 


(10) 

(11) 

(12) 

(13) 

(14) 


^  —~pu  dp  ~  dU,  (7)  =»  d  y  =  ~  ^  [d  P  +  H  (15) 

=  ^d(— --^-  + y  -  ycy)  (16) 

Pcj  =  p0D{D  -  Ucj )  =  pcjUCj{D  -  Ucj)  (17) 

(16)=»i  =  iD[i-(D-U)]  (18) 

The  final  form  is  as  follows 


P 


+  C 


pa-(a,+1) 


ecj  6 


+ 


wp0 


iter-en 


tp 

2 


=  ii  +  ia_r  +  (( 


Pc) 


c; 


—  e 
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(D  -  u)] 


(19) 

(20) 

(21) 

(22) 

(23) 


In  our  implementation,  the  set  of  equations  is  solved  for  a  given  area  expansion,  (r/r0)2 
using  Brent's  method  (ref.  8). 
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ANALYTIC  CYLINDER  MODEL 


Building  upon  Taylor's  long  bomb  model,  one  method  of  including  radial  flow  effects  is  to 
assume  spherical  surfaces  of  constant  thermodynamic  properties  and  mass  flow  in  the  detona¬ 
tion  products  (ref.  9).  The  detonation  products  mass  flow  is  assumed  to  be  in  a  perpendicular 
direction  to  the  spherical  surfaces.  A  diagram  of  a  products  constant  spherical  surfaces  cylinder 
expansion  due  to  high  explosive  detonation  is  presented  in  figure  2.  Again,  it  should  be  noted 
that  flow  velocities  are  relative  to  the  detonation  velocity,  D.  If  constant  detonation  product  prop¬ 
erties  are  assumed  across  spherical  surfaces,  the  following  model  results. 

Mass: 


PcjU  cjAq  pUA 

Axial  momentum: 

Pc;r02  -  Pr 2  =  ^ D2cosQ  -  ^D2  +  pU2r2  -  pCjU2jP o 
Energy: 

pCjt/c;A0  (y-  +  eCj j  +  PCjUcjA0  =  pUA  +  e )  4-  PUA 

Principle  isentrope: 

(<o+i) 


RlP 0  /  V 

p  =  ^ »<«•'  +c(a) 

Taylor  angle: 


1  /I  \ 

’de  =  "Pd(p) 


v  =  2Dsin- 
2 


Spherical  area: 
A  =  7rr2 

sin2  0 


To  achieve  an  easier  computational  form,  the  following  reduction  is  made 

2  (— Y 

(28)  =*  cos0  =  1  —  2  sin2  —  =  1  —  2  (— )  =  1  — 

2  \2DJ  2 


( V/D )2 


(29),  (30)  =>  A  =  2nr 2 — — l 


l- 


tj£L 


nr“ 


(24),  (31) .  p = a£=i  [i  -  fe)2]  -  ag[i  -  fe)2] 


(24) 


(25) 


(26) 


(27) 


(28) 


(29) 


(30) 

(31) 

(32) 
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(33) 


(25)  =>  PcjTq  -  Pr2  =  -=D2 -  Pcluyt }  +  pu2r2 

=K32-^+',©2(,2-^/(,s' 


c 

mpo 


(26)  =>-^  +  ecy+-^  =  ^+  e  +  - 

V  '  2  c>  Pcj  2  p 


=>  —  =  ^+ecy+^-e-- 

2  2  C'  pc;  p 

The  final  form  is  as  follows 

-R iPo 


.-(w  +  l) 


(27 )=>P  =  £iV  p  +  C(^) 

(27)  -  .«>  - .  -  z.^;  ^[(sr  -  (?y 

<36>=*T  =  f  +  S“?  +  e"‘e 

(34)  -  f  =  [P  Q2  -  Pcl  +  p  ©*  U2  -  Pfc;I/», 


c 

mPo 


(34) 

(35) 

(36) 

(37) 

(38) 

(39) 

(40) 

(41) 
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Again,  in  our  implementation,  the  set  of  equations  is  solved  for  a  given  area  expansion, 
(r/r0)2,  using  Brent's  method.  The  Taylor  long  bomb  model  is  appropriate  for  large  charge  to 
mass  ratio  configurations,  where  the  radial  momentum  associated  with  the  detonation  products 
is  relatively  small.  However,  the  spherical  surface  approach  has  been  shown  to  be  much  more 
accurate  for  smaller  mass  to  charge  ratios  without  any  loss  of  agreement  at  larger  mass  to 
charge  ratios  (ref.  9).  It  should  be  recognized  that  this  analytic  modeling  approach  neglects 
initial  acceleration  due  to  shock  processes  (ref.  10)  and  is  therefore  anticipated  to  be  more 
accurate  as  the  initial  shock  process  is  damped  out. 


CYLINDER  WALL  THINNING 

The  analytic  cylinder  models  (fig.  2),  as  expressed,  do  not  consider  the  fact  that  the  cy¬ 
linders  thin  during  radial  expansion.  One  simple  way  to  account  for  this  wall  thinning  is  to  as¬ 
sume  that  the  wall  cross-sectional  area  remains  constant  and  rand  /represents  the  inside  ra¬ 
dius  and  inside  surface  wall  velocity. 

vout  =vrL  ’  rout  =  r?n  +  r02ouC  -  r02ln  (42) 


HIGH  RATE  CONTINUUM  MODELING  COMPARISON 

ALE-3D  (fig.  3)  high  rate  continuum  modeling  was  compared  to  analytic  cylinder  test 
modeling  using  identical  JWLB  equations  of  state  for  TNT  and  LX-14.  The  JWLB  equations  of 
state  were  parameterized  using  JAGUAR  thermochemical  equation  of  state  modeling  (ref.  11). 

T able  1  presents  the  TNT  and  LX-14  JWLB  parameters.  The  1 .2  in.  outer  diameter,  1  in.  inner 
diameter,  10  in.  long  copper  cylinder  was  modeled  using  the  Johnson-Cook  material  model. 
Figures  4  and  5  present  the  comparison  of  the  analytic  cylinder  test  model  to  the  ALE-3D  model¬ 
ing  for  TNT  and  LX-14,  respectively.  The  analytic  cylinder  model  slightly  under  predicts  the  ve¬ 
locities  at  two  and  three  inside  area  expansions,  but  is  in  very  close  agreement  by  six  and  seven 
inside  area  expansions.  This  is  consistent  with  the  fact  that  this  analytic  modeling  approach 
neglects  initial  acceleration  due  to  shock  processes.  Strong  shock  effects  are  typically  observed 
in  the  two  to  three  volume  expansion  region  and  are  pretty  much  damped  out  by  the  six  volume 
expansions,  where  very  close  agreement  between  the  analytic  model  and  ALE-3D  results  are 
observed. 
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Figure  3 

ALE-3D  cylinder  test  modeling  at  10  ps  intervals 
Table  1 

JWLB  parameters  for  various  explosives 


TNT 

LX- 14 

PAX-30 

PAX-30Si 

PAX-30B 

Density  (g/cm3) 

1.6300 

1.8350 

1.9090 

1.8890 

1.8690 

E0  (Mbar) 

0.0657 

0.1032 

0.1376 

0.1358 

0.1441 

D  (cm/s) 

0.6817 

0.8691 

0.8429* 

0.8095* 

0.8584* 

P  (Mbar) 

0.1930 

0.3418 

0.2464* 

0.2601* 

0.2395* 

A1  (Mbar) 

399.2140 

399.1910 

405.3810 

598.5830 

584  5750 

A2  (Mbar) 

56.2911 

52.1951 

14.8887 

71.6056 

0  0130 

A3  (Mbar) 

0.8986 

1 .59892 

1.49138 

1.3717 

1.7536 

A4  (Mbar) 

0.0092 

0.0249 

0.0076 

0.0045 

0.0126 

R1 

28.0876 

27.4041 

13.2982 

25.2998 

12.5850 

R2 

9.7325 

8.4331 

8.0204 

9.9478 

59.7687 

R3 

2.5309 

2.6293 

2.4942 

2.3865 

2.6542 

R4 

6.9817 

0.7498 

0.3566 

0.3112 

0.4086 

All 

58.2649 

68.6476 

66.6542 

0.0028 

0.9895 

AI2 

6.1981 

6.7497 

5.7776 

13.0069 

2.3314 

BI1 

2.9036 

4.1338 

3.1440 

3091.9200 

0.9997 

BI2 

-3.2455 

-4.4607 

-3.2552 

-10.0596 

-0.8634 

RI1 

25.5601 

26.2448 

25.5996 

12.3214 

21.5439 

RI2 

1.7034 

1.6977 

1.7099 

2.0342 

1.2931 

*  Eigenvalue  weak  point  detonation  state  (not  the  Chapman-Jouguet  state) 
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Wall  Velocity  (cm/us)  Wall  Velocity  (cm/us) 


Figure  4 

TNT  Analytic  model  versus  ALE3D  modeling 


Figure  5 

LX-14  Analytic  model  versus  ALE3D  modeling 
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ALUMINUM  (Al),  SILICON  (Si),  AND  BARON  (B)  BASED  COMBINED  EFFECTS  EXPLOSIVES 


The  combined  effects  aluminized  explosives  PAX-29,  PAX-30,  and  PAX-42  were 
previously  demonstrated  to  achieve  both  excellent  metal  pushing  and  high  blast  energies  in  both 
cylinder  test  and  warhead  configurations.  Modeling  comparisons  using  JAGUAR  predicted 
JWLB  parameters  to  cylinder  test  experiments  were  previously  completed  with  very  good 
agreement.  Figure  6  presents  the  PAX-30  results.  More  recently,  new  semi-metal  combined 
explosive  compositions  based  on  Si  and  B  were  investigated  using  the  JAGUAR  thermo¬ 
chemical  equation  of  state.  As  part  of  this  investigation,  improved  modeling  of  Si,  B,  and  their 
possible  detonation  products  was  completed  using  non-linear  optimization  to  parameterize  with 
available  Hugoniot  data,  diamond  anvil  data,  and  melt  curves  (refs.  12  and  13).  These  new 
semi-metal  based  formulations  show  potential  high  energy  and  high  blast  performance  equal  or 
greater  than  current  combined  effects  explosives  and  also  provide  significant  promise  for 
reduced  sensitivity.  Although  B  has  seen  some  investigation  as  an  additive  for  increased 
performance  (refs.  14  and  15),  Si  has  not  been  investigated  for  increased  performance.  B  has  a 
relatively  high  melt  temperature  (2300°C)  compared  to  Si  (1400°C)  and  Al  (660°C),  which  could 
inhibit  its  ability  to  participate  in  a  detonation  reaction.  However,  the  melt  temperature  of  Si 
drops  considerably  with  pressure  and  is  passivated  with  a  much  thinner  oxide  layer  than 
aluminum.  For  this  reason,  Si  is  actually  predicted  to  have  shorter  characteristic  melt  times  than 
Al  at  a  given  particle  size.  For  this  reason,  Si  is  expected  to  participate  in  a  detonation  reaction 
at  larger  particle  sizes  compared  to  Al,  which  could  provide  reduced  sensitivity.  Eigenvalue 
detonation  theory  was  previously  shown  to  explain  the  observed  behavior  of  the  PAX-29,  PAX- 
30,  and  PAX-42  aluminized  combined  effects  explosives  and  is  indicated  for  Si  and  B  based 
combined  effects  explosives  as  well. 


r-rO 
Figure  6 

PAX-30  cylinder  test  experimental  results  comparison  to  modeling 
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EIGENVALUE  ANALYTIC  CYLINDER  MODEL 


Eigenvalue  detonation  theory  is  observed  for  aluminized  combined  effects  explosives 
and  is  anticipated  for  Si  and  B  based  combined  effects  explosive.  For  this  reason,  it  was  of 
interest  to  develop  a  modified  analytic  cylinder  test  model  that  would  provide  a  description  of  the 
detonation  products  isentropic  expansion  from  the  eigenvalue  detonation  weak  point,  rather 
than  from  the  Chapman-Jouguet  state.  The  JWLB  thermodynamic  equation  of  state  is 
essentially  developed  around  a  principle  isentrope  fit.  It  was  found  that  the  most  straightforward 
method  of  implementation  of  an  eigenvalue  detonation  analytic  cylinder  model  was  to  refit  the 
isentrope  associated  with  the  eigenvalue  weak  point  using  the  equation  27  form.  In  this  way, 
equations  24  through  43  remain  correct,  except  that  the  eigenvalue  weak  point  is  used,  rather 
than  the  Chapman-Jouguet  state.  With  this  approach,  it  is  important  to  realize  that  the  weak- 
point  isentrope  fit  is  not  the  same  as  the  principle  isentrope  fit.  The  final  form  is  as  follows 
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Figures  7  through  9  compare  predicted  cylinder  test  results  using  the  eigenvalue  analytic 
cylinder  test  model  for  Al,  Si,  and  B  based  combined  effects  explosives.  Excellent  agreement 
between  the  analytic  cylinder  test  model  and  ALE-3D  cylinder  test  modeling  is  achieved.  The 
compositions  are  HMX  based  with  15%  Al  (PAX-30),  15%  Si,  and  10%  B.  Table  1  presents  the 
JAGUAR  derived  JWLB  parameters  used  in  ALE-3D  for  the  Al  (PAX-30),  Si  (PAX-30Si),  and  B 
(PAX-30B)  based  compositions.  The  JWLB  and  eigenvalue  detonation  weak  point  isentrope 
equation  43  was  parameterized  using  JAGUAR  thermochemical  predicted  behavior.  The  results 
for  these  explosives  show  only  a  very  small  reduction  of  explosive  work  output  for  eigenvalue 
detonations  compared  to  Chapman-Jouguet  detonations.  This  is  due  to  the  fact  that  the 
Chapman-Jouguet  principle  isentrope  and  eigenvalue  weak  point  isentrope  lie  very  close  to 
each  other.  Both  the  Si  and  B  are  predicted  to  have  similar  or  increased  performance  over  LX- 
14,  similar  to  the  Al  based  PAX-30,  which  has  been  experimentally  demonstrated. 
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PAX-30 


Figure  7 

Al  based  combined  effects  explosive  (PAX-30)  analytic  models  versus  ALE-3D 


Figure  8 

Si  based  combined  effects  explosive  analytic  models  versus  ALE-3D 
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PAX-30  30%  Boron 


Figure  9 

B  based  combined  effects  explosive  analytic  models  versus  ALE-3D 


CONCLUSIONS 

A  new  eigenvalue  extended  analytic  cylinder  expansion  model  was  developed  based 
around  isentropic  expansion  from  the  detonation  eigenvalue  weak  point,  rather  than  from  the 
Chapman-Jouguet  state.  Semi-metal  combined  explosive  compositions  based  on  silicone  (Si) 
and  boron  (B)  was  investigated  using  the  JAGUAR  thermochemical  equation  of  state  and  this 
new  eigenvalue  analytic  cylinder  expansion  model.  The  results  for  these  explosives  show  only  a 
very  small  reduction  of  explosive  work  output  for  eigenvalue  detonations  compared  to 
Chapman-Jouguet  detonations.  This  is  due  to  the  fact  that  the  Chapman-Jouguet  principle 
isentrope  and  eigenvalue  weak  point  isentrope  lie  very  close  to  each  other.  These  new  semi¬ 
metal  based  formulations  show  potential  high  energy  and  high  blast  performance  equal  or 
greater  than  current  aluminized  combined  effects  explosives  and  may  also  provide  significant 
promise  for  reduced  sensitivity.  In  particular,  due  to  the  relatively  small  predicted  characteristic 
melt  times,  the  Si  based  explosive  may  provide  high  work  output  at  early  volume  expansions 
using  larger  particle  sizes.  The  use  of  larger  particle  size  may  provide  a  reduction  in  sensitivity. 
In  contrast,  the  high  melt  temperature  of  B  may  inhibit  its  ability  to  participate  in  a  detonation 
reaction,  so  that  predicted  increase  in  early  work  output  may  not  be  achievable. 
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